EQUILIBRIUM FLUCTUATIONS FOR THE 
TOTALLY ASYMMETRIC ZERO RANGE PROCESS 



PATRICIA GONgALVES 

Abstract. We prove a Central Limit Theorem for the empirical measure in the one-dimensional 
Totally Asymmetric Zero-Range Process in the hyperbolic scaling N, starting from the equilib- 
rium measure v p . We also show that when taking the direction of the characteristics, the limit 
density fluctuation field does not evolve in time until iV 4 / 3 , which implies the current across the 
characteristics to vanish in this longer time scale. 



1. Introduction 

In this paper we study the Totally Asymmetric Zero-Range process (TAZRP) in Z. In this 
process, if particles are present at a site x, after a mean one exponential time, one of them jumps 
to x + 1 at rate 1, independently of particles on other sites. This is a Markov process rj. with space 
state N z and configurations are denoted by rj, so that for a site x, rj(x) represents the number of 
particles at that site. For each density of particles p there exists an invariant measure denoted by 
v p , which is translation invariant and such that E Up [rj(0)] = p. 

Since the work of Rezakhanlou in [7] , it is known that for the TAZRP the macroscopic particle 
density profile in the Euler scaling of time, evolves according to the hyperbolic conservation law 
dtp(t,u) + V '<j>(p(t,u)) = 0, where <f>{p) = Since cf> is differentiable, last equation can also 

be written as dtp{t,u) + <fi'(p(t,u))Vp(t,u) — 0. This result is a Law of Large Numbers for 
the empirical measure associated to this process starting from a general set of initial measures 
associated to a profile po, see [7] for details. If one wants to go further and show a Central Limit 
Theorem (C.L.T.) for the empirical measure starting from the equilibrium state v pi one has to 
consider the density fluctuation field as defined below, see (|2.1jl . 

Taking the hyperbolic time scale, the limit density field at time t is just a translation of the 
initial density field. The translation or velocity of the system is given by <j>'(p) — [fjpP which 
is the characteristics speed. If we consider the particle system moving in a reference frame with 
this constant velocity, then the limit field does not evolve in time and one is forced to consider a 
longer time scale. Following the same approach as in [5] we can accomplish the result up to the 
time scale TV 4 / 3 , ie the limit density field does not evolve in time until this time scale. Using this 
approach, the main difficulty in proving the C.L.T. for the empirical measure is the Boltzmann- 
Gibbs Principle, which we can handle by using a multi-scale argument as done for the ASEP in [2], 
but in this case there are some extra computations to overcome the large space state. This result 
implies that the flux of particles through the characteristics speed vanishes in this longer time 
scale. In fact, it was recently proved by pQ that the variance of the current across a characteristic 
is of order t 2 / 3 and this translates by saying that in fact our result should hold till the time scale 
A 3 / 2 . These results should be valid for more general systems than TAZRP or TASEP (see [2j), 
but for systems with one conserved quantity and hyperbolic conservation law. This is a step for 
showing this universality behavior. 

This paper is a natural continuation of [2J and the multi-scale argument seems to be robust 
enough to be able to generalize it to other models and to achieve the conjectured sharp time scale 
A 3 / 2 , this is subject to future work. 
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We remark that all the results presented here, also hold for a more general Zero- Range process, 
namely one could take a Zero-Range dynamics in which the jump rate from x to x + 1 is given by 
g(ri(x)), with g satisfying conditions of definition 3.1 of Chap. 2 of 5J. We could also consider a 
partial asymmetric process, in which a particle jumps from x to x + 1 at rate pg(r)(x)) and from 
x to x — 1 at rate qg{r)(x)), where p + q = 1 and p 1/2 and with g as general as above. The 
results are valid for these more general processes but in order to keep the presentation simple we 
state and prove them for the TAZRP. 

An outline of the article follows. In the second section we introduce the notation and state the 
main results. In the third section we show the Central Limit Theorem for the empirical measure 
in the hyperbolic time scale and the Central Limit Theorem for the current over a fixed bond. 
In the fourth section we use the same approach as in [5] to prove the Central Limit Theorem for 
the empirical measure on a longer time scale and the proof of the Boltzmann-Gibbs Principle is 
postponed to the fifth section. 

2. Statement of Results 
The generator of the one-dimensional TAZRP is given on local functions / : N z -> R by 

^/W = E 1 {n(x)>i}[/(^' x+1 ) ~f(v)], 

x£Z 

where 

7](z), if z 7^ x, x + 1 



V x ' x+1 (z) = { V (x)-l, ifz 



x 



r)(x + l) + l, ifz = x + l 

In order to keep notation the more general as we can, we denote by g(rj(x)) the function l{ r] ( x )>i}, 
which denotes the jump rate of a particle to leave the site x. 

The description of the process is the following. At each site, one can have any integer number 
of particles and for a site x after an exponential time of rate one, one of the particles at that site, 
jumps to the neighboring right site x + 1, at rate 1. Initially, place the particles according to a 
Geometric product measure in N z of parameter p, denoted by v pi which is an invariant measure 
for the process. 

Since the work of Rezakhanlou [7] it is known that taking the TAZRP in the Euler time scaling 
and starting from general initial measures associated to an initial profile po (for details we refer 
the reader to [7]), one gets in the hydrodynamic limit to the hyperbolic conservation law: 



where the flux is given by <fi(p) = 



d t p(t,u)+V<t>(p(t,u)) = 0, 
p 

1+p" 



Fixed a configuration rj, let ir N {rj,du) denote the empirical measure given by 

ir N (rj,du) = ^J2i](x)S^(du) 

x&L 

where 6 U denotes the Dirac measure at u and let ^(77, du) — ir N (r] t , du). 

In order to state the C.L.T. for the empirical measure we need to define a suitable set of test 
functions. For an integer k > 0, denote by Hk the Hilbert space induced by <S(R) (the Schwartz 
space) and the scalar product < /, g >k=< f,K$g >, where < •,• > denotes the inner product 
of L 2 (M.) and Kq is the operator Kq = x 2 — A. Denote by Tt-k the dual of Hk, relatively to the 
inner product of L 2 (R) . 

Fix p and an integer k. Denote by y N the linear functional acting on functions H e 6>(R) as 

y t N (H) = Vn\ < ff.Trf fadu) > -E^ p < H,Tr?(r),du) > 

(2.1) =^£^)(^(*)-^ 
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where < H , tt^ (ij , du) > denotes the integral of a test function H wrt to the measure irf (rj,du). 
Throughout the article the functional above is mentioned as the density fluctuation field of the 
process. Denote by D(R + ,H.-k) (resp. C(R + ,H-k)) the space of f/_fc-valued functions, right 
continuous with left limits (resp. continuous), endowed with the uniform weak topology, by Qn 
the probability measure on D(R + ,H-k) induced by y N and v p . Consider = P Up the p.m. on 
Z)(R + , {0, 1} Z ) induced by v p and r\t speeded up by N and denote by E„ the expectation with 
respect to P v . 

Theorem 2.1. Fix an integer k > 2. Denote by Q be the probability measure on C(R + ,7i_fc) 
corresponding to a stationary Gaussian process with mean and covariance given by 

E Q [y f (H)y s (G)} = X (p) J H(u + <P'( P )(t - s))G(u)du 

for every < s < t and H, G in TCk- Here x{p) = V<ir(r)(0),i/ p ). Then, (Qn)n>i converges 
weakly to Q. 

Last result holds for the TAZRP evolving in any Z d and for the other more general processes 
as pointed out in the introduction. 

An easy consequence of last result is the derivation of the C.L.T. for the current over a fixed 
bond, see [3]. For a site x, denote by J x v x+1 (tN) the current through the bond [x,x + 1], defined 
as the total number of jumps from x to x + 1 during the time interval [0, tN]. 



Theorem 2.2. Fix x e Z, t > and let 

Z t N = -j={j» x+1 (tN) -E Up [j£ x+1 (tN)]}. 



/N 

Then, under P y , 



B 



weakly, where B t denotes the standard Brownian motion. 

In the hyperbolic scaling, the limit density fluctuation field at time t is a translation of the 
initial one and removing from the system this translation velocity, it does not evolve in time and 
one is forced to go beyond the hydrodynamic time scale. With this in mind, let r\. be evolving 
in the time scale N 1+1 ', with 7 > 0, fix p and remove the translation velocity by redefining the 
density fluctuation field on H £ <S(K) by: 



As above, let Q^ N be the probability measure on D{R+,H- k ) induced by y N -~< and v p , let 
pjv,7 _ ^ c p m Qn £j^+ ; |q j i} z ) induced by v p and 77* speeded up by N 1+ ~> and denote 
by Wl p expectation with respect to PJ . Then 

Theorem 2.3. Fix an integer k > 1 and^ < 1/3. LetQ be the probability measure on C(R. + ,7i 
corresponding to a stationary Gaussian process with mean and covariance given by 



E Q [y t (H)y s (G)} = x (p) j H(u)G(u)du 

for every s,t > and H, G in Tlk- Then, (Qjy)jv>i converges weakly to Q. 

The main problem to overcome when showing last result is the Boltzmann-Gibbs Principle, 
which we can prove for 7 < 1/3 using a multi-scale argument as for the ASEP in [2]. 

Theorem 2.4. (Boltzmann-Gibbs Principle) 
Fix 7 < 1/3. For every t > and H £ S(R), 

2 



lim EZ 



= 0, 



4 
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where 

V g (v(x)) = g(r](x)) - cf>{p) - cf>'(p)[v(x) - p], 

and 4>(p) = EvMrtO))]. 

Now we define the current of particles across a characteristic. Let J^i n (tN) be the current 
through the bond [vf ,vf + 1] (where vf = x + [</>' (p)tN 1+ ~<]) defined as the number of particles 
that jump from vf to vf + 1, from time to tN 1+1 : 



4 n (tN) = £ (rH(y + vf) - Vo (y + xj) . 



y>l 

As a consequence of last result, it holds that: 

Proposition 2.5. Fix t > 0, a site x G Z and 7 < 1/3. Then, 



lim 



N 



= 0. 



3. Density Fluctuations for the Hyperbolic Scaling 



3.1. Equilibrium Fluctuations. Fix a positive integer k, denote by 21 the operator 0'(p)V 
defined on a domain of L 2 (R) and by {T t ,t > 0} its semigroup. The theorem follows as long as 
we show that {Qn)n>i is tight and characterize the limiting measure Q. 
Fix H £ 5(E), then 

M t N ' H = 3^(2Z) - y»(n) V*ir(£) ff fo.(*))<fa 
is a martingale with respect to the natural filtration with quadratic variation given by 



ij 2 (v"ff (^)) 2 [<?(%(*)) + g{ris{x + 1)) 



where \7 N H denotes the discrete derivative of i/. The integral part of the martingale can be 
written as 



ds. 



by using the fact that J2 X £Z^ = 0- The following result allows to replace g(n s (x)) — 4>{p) 

by 4>' {p)[rj s {x) — p] and allows to recover the density fluctuation field inside the integral part of 
the martingale. 

Theorem 3.1. (Boltzmann-Gibbs Principle) 
For every H £ 5(R) and every t > 0, 



lim E„ 

N-*oa " 



1 



i? 



(|f)^M*)K' 



= 0. 



The proof of last result follows the same lines as for the Symmetric Zero-Range Process in [5] 
and for that reason we have omitted it. For the same reason the following results are just stated 
but their proofs follow the same lines as for the ASEP in [2]: (Qn)n>i is a tight sequence, the 
limiting measure Q is supported on fields y. such that for a fixed time t and a test function H, 
yt(H) = yo(TtH) where TtH(u) = H (u + tf>' (p)t) and 3^0 is a Gaussian field with covariance given 
by E Q (y (G)y (H)) = x (p) <G 1 H>. 
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3.2. Central Limit Theorem for the Current over a fixed bond. Now we give a sketch 
of the proof of Theorem 12.21 in which we need to show the convergence of finite dimensional 
distributions of Zf j \Jx{p)4>'{p) to those of Brownian motion together with tightness. 

We start by the convergence of finite dimensional distributions, namely, we show that for every 
k > 1 and every < t\ < t 2 < ■■ < tk, (Z^ , .., Z^) converges in law to a Gaussian vector 
(Z tl , .., Z tk ) with mean zero and covariance given by EQ[Z t Z s ] = x(p)4>'{p) s provided s < t. 

Recall that J^ 1Q (tN) is defined as the total number of jumps from the site —1 to during the 
time interval [0,iiV]. Since 

J N h0 (tN) = J2(vt(x)-Vo(x)] 



x>0 



the current can be written in terms of the density fluctuation field evaluated on Hq, the Heaviside 
function H (u) = l[ 0)0o )(«): 

±={j%(tN) - E Up [J» h0 (tN)}} = y»(H ) - y»(H ) t 



<N 

Approximating Hq by {G n ) n >i such that G n (u) = (1 — — ) + lro )00 )( u ); then 
Proposition 3.2. For every t > 0, 



lim E„ 

n — *+oo 



rJV,ii,Vi - (y?(G n ) - y»(G n )) 



N 



2 



uniformly in N . 



The convergence of finite dimensional distributions is an easy consequence of last result together 
with Theorem 12.11 see [3]. 

Now, it remains to prove that the distributions of Zf / y/x{p)4>' (p) are tight. For that, we 
can use the same argument as in Theorem 2.3 of [2] that relies on the use of Theorem 2.1 of 
[5] with the definition of weakly positive associated increments given in [S]. One can follow the 
same arguments as those of Theorem 2 of [4] to show that J-i,o(i) has weakly positiive associated 
increments with the definition in [9] , see [2] . In order to conclude the proof it remains to note that 

lim jE Up [J- ll0 (t)] = cr 2 , 

which follows by Theorem 3 of 0]. 



4. Density Fluctuations for a longer time scale 



Fix a positive integer k and let H{u) — H(u — <j)' (p)tN J ). Recall the definition of (QX-)jv>i 
and note that following the same computations as in [5] it is easy to show that the sequence is 
tight. Now we compute the limit field, by fixing H € <S(R) such that 

Jo V x£Z 



is a martingale and whose quadratic variation is given by 



7 N TT N 



g(Va(x)) +g(vs(x + i)) 



ds . 



If 7 < 1, Mt' H vanishes in L' 2 {P2 ) as N 



Using the Botzmann-Gibbs Principle, whose 



proof is sketched in the next section, the integral part of the martingale M t ' vanishes in L 2 (P^ ) 
as N — > +oo which in turn implies that if Q is one limiting point of (Qn)n, the limit density 
fluctuation field satisfies yt(H) = yo(H), where 3^0 is a Gaussian field with covariance given by 
E Q (y (G)y (H)) = x(p) <G,H>. 



(i 
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5. BOLTZMANN-GlBBS PRINCIPLE 

In this section we prove Theorem l2.4l Since we are going to follow the same steps as in Theorem 
(2.6) of [2] we just remark the fundamental differences between the proofs. 

To start fix an integer K and a test function H £ <S(K). We divide Z in non overlapping 

intervals of length K, denoted by {Ij,j > 1} and by summing and subtracting Hy^j, where yj 

is some point of Ij , we can bound the expectation appearing in the statement of the Theorem by 



2E7, 



ATT 



+2E? 



££[ 

t _/v 7 



" X V j>lxElj 



V g (ri s (x))ds 



£#(!)£ F ^ x )) rfs 



j>l x£lj 

The first expectation is easily handled, since by Schwarz inequality and the invariance of v p it 

can be bounded by Ct 2 N 2 ^\ \H'\ || (f ) and vanishes as long as KN 1 1 — > when N — ► +oo. 

In order to treat the remaining expectation we bound it from above by 
(5.1) 



2m 

Vp 

where 



j>l v j>l xG/j 



M, 



vuAv) = E ^w*)) - E ( E w*)) 

and = cr^xe/, vi*))- 

Lemma 5.1. For every H £ S(M) and every t > 0, if K 2 N 1 ^ 1 — > as N — > +oo, i/ien 



lim 



= 0. 



Proof. By Proposition Al.6.1 of [5] and by the variational formula for the _ff_!-norm the expec- 
tation above is bounded by 

CtY, sup {2 / ELh (|Vi,i^)%K(*?) - < h,-£f.h > p }, 

where £ is the Symmetric dynamics restricted to the set Ij, namely: 

= E ^hwyil/M-ZW]' 



where 



[as— 1/|=1 



77(2:), iHz^x,y 
rf' v (z) = { rj(x) — 1, if 2 = a; 

77(2/) + 1, if z = y 



For each j and Aj a positive constant, it holds that 



and taking Aj 



H(w)\) > t nc whole expectation becomes bounded by 

«E^ 2 (|)<^>(-4)- 1 ^ 
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By the spectral gap inequality for the Symmetric Zero- Range process (see [B]) last expression can 
be bounded by 

N 2 \N, 



The proof of the Lemma ends if we show that Var(V\ j j3 , v p ) < KC, since it implies that the 
expectation in the statement of the lemma to be bounded by Ct^-(K + l) 2 ||-ff||2 and vanishes 
as long as K 2 N 1_1 -> when N -> +00. □ 

Remark 5.1. Here we show that Var(Vi t j ig ,u p ) < KC. Since Var(Vij, g ,iy p ) < E Up [Vi j, s ] 2 and 
by the definition of Vij ig we have that 

Var(V ldfg ,u p ) <^, P [E V Mx))~Eu p lJ2 V M*))\ M i\ 2 



xeij 



xeij 



By the definition of V g (rj) last expression can be written as 

Eu„ [ £ (a(v(x)) - Hp) ct>'( P )[ v (x) - P ])-J2 Mp) *Hp) E <I>'<J>M p\ 

xeij x ^Ij x ^Ij 

where <fcj(p) = E Vp [g(rj)\Mj\ and rff — ^ ^ l£j v( x )- On the other hand, by summing and sub- 
tracting Hvf) — E v x [g(r))], where v^k is the Bernoulli measure with density rfc , last expression 
can be bounded by 

2E - P [ E {aW*)) - Hp))] 2 + 2£, p [ £ <f/(j>)[r,(x) - p] 



xeij 



xeij 



2E », [ E (Mp) - Hvf))] 2 + 2£, p [ £ (Hvf) Hp) - H{ P )[vf p\ 



xeij 



xeij 



Now we treat each expectation separately. 

For the first and the second one, since {rj{x)) x are independent under v p , it is easy to show that 

E »„ [ E (9(V( X )) - HP)) ] 2 < CKVar(g, v p ) 

xeij 

and 

E„ p [Y,H(p)lv(x)-p}] 2 < CKVar(r,(0),u p ). 

xeij 

On the other hand, to treat the third expectation one can use the equivalence of ensembles (see 
Corollary A2.1.7 of [S]J which guarantees that \(f>j(p) — Hvf)\ — ^~k^~ while for the last one, one 
can use Taylor expansion to have 



Hvf)-Hp)-H(p)(vf-p) 



vf-P 



= 0(K- 2 ). 



Putting these arguments all together one gets to the bound KC. 



To conclude the proof it remains to bound the expectation on the right hand side of (|5.ip ■ For 
that, fix an integer L and take disjoint intervals of length M — LK, denoted by {Ii,l > 1} and 
write it as: 



EI 



AH 



E£*($MEW*)) 



'>* J67, 



xeij 



Mj ) ds 



By summing and subtracting H ( ^ ) , where zi denotes one point of the interval // , last expectation 
can be bounded by 



* ATT 
N 



EE[ 

i>i 3 eh 



H 



H 



xeij 
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Following the same arguments as above it is easy to show that the first expectation vanishes if 
L 2 KN 2l ~ 2 — ► as N — > +00. For the second one, sum and subtract E^J^xeii Vgi'ti 30 )) M 
where M; = o~( ^2 x eh 7 l( x )j an d bound it by 



2E2 
where 



' ^ 1 ' ^ 1 X £ 



M; Ids 



w?) = E ( E w^K) - E ( E w*)) 



Mi 



xEli 



Lemma 5.2. For every H S 5(R) cmd every i > 0, if L 2 KN" f 1 — > as — > +oo, i/ien 



lim E7 



N 1 



E^(|)w%K< 
01 



0. 



Proof. Following the proof of Lemma (|5.1[) , the expectation becomes bounded by 

Ci^ sup {2 f ^ H ^)v 2 ,iAv)h(vH(dv)-N 1+ ~< <h 1 -£ s i h> p }. 

Using an appropriate ^4/ and the spectral gap inequality, last expression is bounded by 

Ct E ^ (f ) ( M + ^VariV^g, v P ). 
i>i 

Now, the proof ends as long as Var(V2,i t g, v p ) < LC, which is proved below. 



□ 



Remark 5.2. Here we show that Var{V2,i t g,v p ) < LC . Since Var(V2,i, g , v p ) < E Vp \V2,i,g] 2 and 
by the definition of V2,i,g we have that 

2 



Var{V 2 ,i, g ,v p ) <^ P [E S (E VM^Mj) - ^( E V M*)) 



jeh 



Mi 



xEh 



By the definition of V g (rj) and the notation introduced above, one can write last expression as 

E Up [ (K4>M ~ Kcb(p) - <j>'{p)K[ n f - p\) - MMp) ~ Mcb(p) - Mtf/{p)[rf - p] 
jeh 

where (j>i(p) — E Up \g{rf)\Mi\ and rjf' 1 = jj ^2 xeIl {rj{x) — p). Last expression can be written as 
1 _ / n -m2 



M { m E ( K ^(p) ~ K ^p) - #(p)K[tf - pi) - mp) - m - ^'{pU 1 P}} 



J eh 




= E u 


M 


M 2 




= —E„_ 


L 


- p 



jeh 



jeh 



By the independence of the random variables (rj(x)) x under v p and the Central Limit Theorem, 
last expectation is of order 

2 



E„ 



4>M-<t>(p)-V{p)K -p\ 



which we can bound by 



2E„ 



<Pj (p) - 



2E V 



<i>(v?)-M-^(p)K -p] 
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By the equivalence of ensembles the expectation on the left hand side is bounded by K 2 . For the 
other, use Taylor expansion to have 



E„ 



<P(vf) ~ m ~ 4>'{p)\rif -P] ~ E Up [r,f - pf = 0(K- 2 ) 



This finishes the proof of the remark. 



The proof of Boltzmann-Gibbs Principle 

Following the same arguments as before, take n sufficiently big for which in the n-th step of the 
proof we have intervals, denoted by {Ip,p > 1 >} of length K n = TV 1-7 . At this stage it remains 
to bound: 



p>i xei" 



M™)ds 



where for each p, z p is one point of the interval 7™ and M™ = cr(^2 xe in r](x) 
Since v p is an invariant product measure, last expectation can be bounded by 



(5.2) 



p>i 



p 



Remark 5.3. Here we show that Eu p (E U/ ,lYlxei n ^(^C 2 -)) 
By the definition ofV g , the expectation above is equal to 



Ml 



= 0(1). 



E " P (X ( E (p) - <t>(p) - MM" - p] 



and bounded from above by 

2E " P (X ( E (P) - "))) 2 + 2 ^p (X ( E (*fan ") - <I>{P) ~ <P'{pM n - P] 



xel" 



where 4>K n (p) = E Vp [g(rj)\Mp] and rj^ n = Yl x el n 7 ?( :E )- Now the result follows if one applies 
equivalence of ensembles to the expectation on the left hand side and Taylor expansion to the 
expectation on the right hand side. 

This implies (|5.2|) to be bounded by which vanishes as N — > +oo since 7 < 1/3. 

Remark 5.4. Here we give an application of the Boltzmann-Gibbs Principle for a linear functional 
associated to the one- dimensional Symmetric Zero-Range process, in the diffusive scaling. Consider 
a Markov process rjtN 2 with generator given by 



x,y£I* 
\x—y\—\ 



with n x ' v as defined in the proof of Lemma \5.1\ If one repeats the same steps as done in the proof 
of Theorem \2.J\ it is easy to show that: 



= 0. 



Corollary 5.3. Fix t > and [3 < 1/2, then 



lim Ej, 

N^oo 



So, in order to observe fluctuations for this field one has to take (3 > 1/2. 



in 
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5.1. Current through the characteristics speed. As in the hyperbolic scaling, Proposition 
(|2.5p is a consequence of: 

Proposition 5.4. For every t >Q and 7 < 1/3: 



uniformly over N . 

The proof of this result follows the same lines as the proof of Proposition 9.4 in [2] and for that 
reason we have omitted it. 
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